Electron localization of linear symmetric molecular ion H 3 2 +
Jia Zheng-Mao, Zeng Zhi-Nan , Li Ru-Xin
State Key Laboratory of High Field Laser Physics, Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Shanghai 201800, China

 

† Corresponding author. E-mail: zhinan_zeng@mail.siom.ac.cn

Abstract

Electron localization in the dissociation of the symmetric linear molecular ion is investigated. The numerical simulation shows that the electron localization distribution is dependent on the central frequency and peak electric field amplitude of the external ultrashort ultraviolet laser pulse. When the electrons of the ground state are excited onto the 2p by a one-photon process, most electrons of the dissociation states are localized at the protons on both sides symmetrically. Almost no electron is stabilized at the middle proton due to the odd symmetry of the wave function. With the increase of the frequency of the external ultraviolet laser pulse, the electron localization ratio of the middle proton increases, for more electrons of the ground state are excited onto the higher 3p state. 50.9% electrons of all the dissociation events can be captured by the middle Coulomb potential well through optimizing the central frequency and peak electric field amplitude of the ultraviolet laser pulse. Besides, a direct current (DC) electric field can be utilized to control the electron motions of the dissociation states after the excitation of an ultraviolet laser pulse, and 68.8% electrons of the dissociation states can be controlled into the middle proton.

1. Introduction

Coherent controls of electrons and fragments in chemical reactions and photoelectron processes have attracted a great deal of interest. One of the main goals is to find a way to selectively break and form molecular bonds in photochemical reactions.[13] For searching for the underlying mechanism in realizing electron localization control in the dissociation process, several solution routes have been proposed, which include the mixture of the 1s and 2p states,[4] the interference between the 2p and states,[5] the superposition of the 2p and 2p states,[6] etc. In addition, the other quantum coupled equations[79] or laser-induced Stark shift effect[1012] have also been used to reveal the dissociation control mechanism.

As is well known, the linear molecule in does not exist in the field-free case. However, the one-electron linear H can be stabilized at high intensities and frequencies due to high nonlinear electron-field interaction.[13] The linear molecular ion can also exist in a strong magnetic field, and the direction of the molecular axis is parallel to the direction of the magnetic field.[14] Furthermore, the collision between and or the dissociative ionization of can obtain too.[15]

Compared with the simplest double-well molecular ion and its isotopes, the linear molecular ion includes three nuclei along the molecular axis. The main goal of this paper is to enhance the electron localization ratio of the dissociation states of the middle proton. The simulation results show that the electron localization ratio of the middle proton is dependent on the central frequency and peak electric field amplitude of the external linearly polarized ultrashort ultraviolet (UV) laser pulse. When the central frequency of the UV laser pulse is 0.2 a.u., about 31.8% electrons of the ground 1s state are excited into the exciting 2p state through a one-photon process. There exists a symmetric electron localization distribution, as a result of the symmetric distribution of the Coulomb potential wells of the linear molecular ion . The most electrons of the dissociation states are stabilized at the protons on both sides, and almost no electron is localized at the middle proton due to the odd symmetry of the wave function of the 2p state. With the increase of the central frequency of the UV laser pulse, more electrons of the ground state are excited onto the higher 3p state by a three-photon process. The electron localization ratio of the dissociation states of the middle proton increases from 0.3% to 50.9%, by optimizing the central frequency and peak electric field amplitude of the external UV laser pulse. Besides, a direct current (DC) electric field can be used to steer the electron motion after the excitation of a UV laser pulse. The symmetric electron localization distribution is broken seriously, as a result of the dressing effect of the DC electric field. The electrons of the dissociation states of the dressed-down potential well move in the direction opposite to the dc electric field force, and are captured by the middle potential well. 68.8% electrons of all the dissociation states can be steered onto the middle proton with the variation of the amplitude of the DC electric field.

2. Simulation model

A reduced-dimensional model for the symmetric linear molecular ion is utilized in the calculation. The molecular axis is assumed to be parallel to the polarization direction of the UV laser field. The middle proton P1 is set to be in the center of the coordinate system. The inter-nuclear distances from P3 and P2 to P1 are . In other words, the UV laser pulse is assumed to have no effect on the motions of the protons. Here P2 and P3 are the protons located on each side of P1. The R is the relative inter-nuclear distance between P2 and P3, and more details can be found in Fig. 1. Then we can adopt the one-dimensional non-Born–Oppenheimer time-dependent Schrödinger equation (TDSE) to conduct the simulation. The corresponding TDSE can be written as ( in atomic unit (a.u.), which are used throughout the paper unless otherwise stated)[1618]

(1)
where is the field-free Hamiltonian of the system, refers to the soft-core Coulomb interaction, indicates the interaction of the particle with the external laser pulse and z represents the electronic coordinate with respect to the center of the coordinate system.

Fig. 1. Geometrical setting for the symmetric linear molecular ion along the z axis. The middle proton P1 is stabilized in the center of the coordinates. The inter-nuclear distances from P3 and P2 to P1 are , where R is the relative inter-nuclear distance between the protons of P2 and P3.

For our model, the kinetic energy in Eq. (1) reads

(2)
then the Coulomb potential of the system can be expressed as[19]
(3)
Here and are the electron mass and proton mass ( , , respectively. The soft-core parameters here are , and .

The interaction of the system with the external laser field can be described as

(4)
Here refers to the electric component of the UV laser field.

The time dependent electric field component of the UV pulse is defined as . Here E0 is the peak electric field amplitude in atomic units, T is the pulse duration, and ω is the central frequency of the UV pulse. The three dissociation channels (see Fig. 2), i.e., the electron localization probabilities of the left, middle, and right protons are defined, respectively, as

(5)
Here corresponds to the boundary range of the R axis, and is the final wave function of the system. In the simulation, is taken as 92.6 fs after the on-set of the UV pulse, when the probabilities with which the electron are localized, respectively, on the three protons (left , middle , and right ), are all stable. The electron localization ratios of the dissociation states of these three protons are set to be , respectively.

Fig. 2. Configuration space regions associated with different dissociation channels. The four dotted curves marked with I, II, III, and IV can be expressed as , , , and , respectively.
3. Simulation results and discussion

In the simulation, a UV pulse with an intensity of W/cm2 and a pulse duration of 10.6 fs is used to excite the electron wave packet onto the dissociative states. Figure 3(a) shows the plots of electron localization probabilities versus central frequency ω of the UV laser pulse for different localization ratios of the three protons. From this figure one can find that there exists a symmetric electron localization distribution with due to the symmetric distribution of the Coulomb potential well of the linear molecular ion . The electron localization probabilities and localization ratios of these three protons are all dependent on ω. When the central frequency is a.u., the central wavelength of the UV laser pulse is 228 nm. The electron localization probabilities of these three protons are , and , respectively. Most electrons of the dissociation states are stabilized at the protons on both sides as shown in Fig. 4(a). When the central frequency of the UV laser pulse is 0.2 a.u., about 31.8% electrons of the ground 1s state are excited onto the exciting 2p state through a one-photon process, see Fig. 3(b).[20] For the electrons of the 2p state, most of them are stabilized at the protons on both sides, and almost no electron is localized at the middle proton due to the odd symmetry of the wave function, which can be seen from Fig. 4(c). At the end of the simulation, the electron localization ratios of the dissociation states of these three protons are , and , respectively. Only 0.3% electrons of the dissociation states are localized at the middle proton.

Fig. 3. (color online) (a) Variations of electron localization probabilities with central frequency ω of the UV laser pulse for different localization ratios of the three protons, where the UV laser pulse intensity is W/cm2. (b) Variations of lowest four electronic energy levels of field-free linear symmetric molecular ion with R.
Fig. 4. (color online) Upper panels: snapshots of the common logarithm of the electron-nuclear probability density distribution taken at the end of the dissociation, with the central frequencies of the UV laser pulse being (a) 0.2 a.u. and (b) 0.292 a.u. Lower panels: the wave functions corresponding to (c) 2p and (d) 3p states, respectively.

With the increase of ω, more electrons of the dissociation states are captured by the middle proton. For example, when the central frequency is a.u., the electron localization probabilities of these three protons are , and , respectively, which can be seen from Fig. 4(b). The 47.2% electrons of the dissociation states are stabilized at the middle proton at the end of the simulation. When the central frequency of the UV laser pulse is a.u., some electrons of the ground state 1s state are excited onto the higher 3p state through a three-photon process (see Fig. 3(b)). Electrons of the 3p state can be captured by the middle proton, because the wave function includes two peaks near z=0.0a.u. as shown in Fig. 4(d).

When the central frequency of the UV laser pulse is a.u., the electron localization ratio of the middle proton decreases due to the decrease of the electron localization probability of the 3p state.

From Figs. 4(a) and 4(b) one can find that the electron localization starts from about R=7.2 a.u. Thus we use a.u. in the definition of , see Eq. (5).

Figure 5 depicts the plots of electron localization probabilities versus peak electric field amplitude E0 of the UV laser pulse for different localization ratios of the three protons. The central frequency of the UV laser pulse is 0.292 a.u.. When E0 is low, most electrons of the dissociation states are localized at the protons on both sides. For example, when the peak electric field amplitude of the UV laser pulse is a.u., , and . The electron localization ratios of the left and right protons are , only 0.9% electrons of the dissociation states are stabilized at the middle proton which is because when the peak intensity of the UV pulse is low, the one-photon process is dominant. The electrons of the ground state can be excited only onto the exciting 2p state through a one-photon process. For the electrons of the 2p state, most of them are stabilized at the protons on both sides, and almost no electron is localized at the middle proton due to the odd symmetry of the wave function. With the increase of E0, more electrons are localized at the middle proton, which is because when the central frequency of the UV laser pulse is 0.292 a.u., the electrons of the ground state can be excited onto the higher 3p state through a three-photon process. The electrons of the 3p state can be captured by the middle potential well. The excitation ratio of the three-photon process is dependent on the peak electric field amplitude E0 of the UV laser pulse seriously. For example, when the peak electric field amplitude is 0.0475 a.u., the electron localization probabilities of these three protons are , and , respectively. The 50.9% electrons of all the dissociation events are stabilized at the middle proton.

Fig. 5. (color online) Plots of electron localization probabilities versus peak electric field amplitude E0 of the UV laser pulse for different localization ratios of the three protons. The central frequency ω of the UV laser pulse is 0.292 a.u.

With the further increase of E0, the electron localization probabilities of these three protons all drop. More electrons are ionized and escape away.

The other way to enhance the electron localization ratio of the dissociation states of the middle proton is to use a DC electric field, whose polarization is assumed to be parallel to the molecular axis, to steer the electron motion after the excitation of an ultrashort UV laser pulse. For the UV pulse, the pulse duration, central laser frequency, and intensity are 10.6 fs, 0.2 a.u., and W/cm2, respectively. Figure 6(a) shows the variations of the electron localization probabilities with the amplitude of the dc electric field for different electron localization ratios of these three protons. When approaches 0.0, there appears a symmetric electron distribution with , and almost no electrons of the dissociation states are stabilized at the middle proton due to the symmetric distribution of the Coulomb potential wells of the linear symmetric molecular ion and the odd symmetry of the wave function of the 2p state.

Fig. 6. (color online) (a) Plots of electron localization probabilities versus the amplitude of DC electric field for different localization ratios of the three protons, in which a DC electric field is utilized to steer the electron motions of the dissociation states after being excited by a UV laser pulse with a central frequency and intensity of 0.2 a.u. and W/cm2, respectively. (b) Snapshot of the common logarithm of the electron-nuclear probability density distribution taken at the end of the dissociation, with the amplitude of the dc electric field being a.u.

With the increase of the amplitude of the dc electric field , the symmetric electron localization distribution is broken seriously, as a result of the dressing effect of the DC electric field. The ionized electrons escape away along the DC electric field force, see Fig. 6(b), while most of the electrons in the dissociation states, which are localized at the left proton subjected to a single UV laser pulse, move in the direction opposite to the DC electric field force, and are stabilized at the middle proton.[21] When a.u., the electrons on the right proton start to be ionized away while the electrons of the middle proton are stabilized due to a deeper potential well. When the amplitude of the DC electric field is a.u., the electron localization probabilities of these three protons are , , and , respectively, which can be seen from Fig. 6(b). The 68.8% electrons of the dissociation states are stabilized at the middle proton at the end of the simulation.

With the further increase of the amplitude of the dc electric field, the electron localization ratio of the middle proton decreases, for a further dissociation is induced.

4. Conclusions

In this work, the electron localization in the dissociation of the symmetric linear molecular ion is investigated. The numerical simulation shows that the electron localization ratio of the middle proton is dependent on the central frequency and peak electric field amplitude of the external UV laser pulse. When the electrons of the ground state are excited onto the 2p by a one-photon process, almost no electron is stabilized at the middle proton due to the odd symmetry of the wave function. With the increase of the central frequency of the UV laser pulse, more electrons of the dissociation states are stabilized at the middle proton, because more electrons of the ground state are excited onto the higher 3p state through a three-photon process. The electron localization ratio of the middle proton can be raised to 50.9% by optimizing the central frequency and peak electric field amplitude of the UV pulse. Besides, a DC electric field can be utilized to control the electron motions of the dissociation states after being excited by an ultrashort UV laser pulse. The symmetric electron localization distribution is broken seriously due to the dressing effect of the DC electric field. The electrons of the dissociation states, which are localized at the potential well of one side, move in the direction opposite to the DC electric field force and are stabilized at the middle proton. As the amplitude of the DC electric field changes, 68.8% electrons of all the dissociative events can be controlled onto the middle proton.

Reference
[1] Gaire B McKenna J Zohrabi M Carnes K D Esry B D Ben Itzhak I 2000 Phys. Rev. A 85 023419
[2] Roudnev V Esry B D 2007 Phys. Rev. A 76 023403
[3] Roudnev V Esry B D Ben Itzhak I 2004 Phys. Rev. Lett. 93 163601
[4] McKenna J Anis F Sayler A M Gaire B Johnson Nora G Parke E Carnes K D Esry B D Ben Itzhak I 2010 Phys. Rev. A 85 023405
[5] Fischer B Kremer M Pfeifer T Feuerstein B Sharma V Thumm U Schröter C D Moshammer R Ullrich J 2010 Phys. Rev. Lett. 105 223001
[6] He F 2012 Phys. Rev. A 86 063415
[7] Kling M F Siedschlag Ch Znakovskaya I Verhoef A J Zherebtsov S Krausz F Lezius M Vrakking M J J 2008 Mole Phys. 106 455
[8] Kling M F Siedschlag Ch Verhoef A J Khan J I Schultze M Uphues Th Ni Y Uiberacker M Drescher M Krausz F Vrakking M J J 2006 Science 312 246
[9] Sansone G Kelkensberg F Pérez Torres J F Morales F Kling M F Siu W Ghafur O Johnsson Swoboda P M Benedetti E Ferrari Lépine F F Sanz Vicario J L Zherebtsov S Znakovskaya I L’Huillier A Ivanov M Yu Nisoli M Martín F Vrakking M J J 2010 Nature 465 763
[10] Dietrich P Ivanov M Yu Ilkov F A Corkum P B 1996 Phys. Rev. Lett. 77 4150
[11] Niikura H Corkum P B Villeneuve D M 2003 Phys. Rev. Lett. 90 203601
[12] Niikura H Villeneuve D M Corkum P B 2004 Phys. Rev. Lett. 92 133002
[13] Zuo T Bandrauk A D 1995 Phys. Rev. A 51 R26
[14] Turbiner A V López Vieyra J C 2006 Phys. Rep. 424 309
[15] Bandrauk A Barmaki S Kamta G L 2007 Phys. Rev. Lett. 98 013001
[16] Jia Z M Zeng Z N Li R X Xu Z Z Deng Y P 2014 Phys. Rev. A 89 023419
[17] Jia Z M Zeng Z N Li R X Xu Z Z 2014 Chin. Phys. B 23 083201
[18] Jia Z M Zeng Z N Li R X Xu Z Z Deng Y P 2015 Chin. Phys. B 24 013204
[19] Bian X B Peng L Y Shi T Y 2008 Phys. Rev. A 78 053408
[20] Herzberg G 1953 Molecular Spectra and Molecular Structure I Spectra of Diamotic Molecules 2 New York National Astronomical Observatories, CAS and IOP Publishing Ltd. 321 324
[21] Jia Z M Zeng Z N Li R X Xu Z Z 2016 J. Phys. B 49 215604